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ABSTRACT 


Recently several second order closure models have been proposed for closing the second 
moment equations, in which the velocity-pressure gradient (and scalar-pressure gradient) 
tensor and the dissipation rate tensor are the two of the most important terms. In the 
literature, these correlation tensors are usually decomposed into a so called rapid term and 
a return-to-isotropy term. Models of these terms have been used in global flow calculations 
together with other modeled terms. However, their individual behavior in different flows 
have not been fully examined because they are un-measurable in the laboratory. Recently, 
the development of direct numerical simulation (DNS) of turbulence has given us the 
opportunity to do this kind of study. With the direct numerical simulation, we may use 
the solution to exactly calculate the values of these correlation terms and then directly 
compare them with the values from their modeled formulations (models). In this paper, we 
make direct comparisons of five representative rapid models and eight return-to-isotropy 
models using the DNS data of forty five homogeneous flows which were done by Rogers et 
al. (1986) and Lee et al. (1985). The purpose of these direct comparisons is to explore the 
performance of these models in different flows and identify the ones which give the best 
performance. The paper also describes the modeling procedure, model constraints, and 
the various evaluated models. The detailed results of the direct comparisons are discussed, 
and a few concluding remarks on turbulence models are given. 


1 


1. Introduction 


Homogeneous turbulent flows have been used for studying many important phenomena of 
turbulence. In these flows, the mean flow field is decoupled from the dynamical equations of 
turbulence and, therefore, the detailed turbulence structure and related turbulent transfer 
have no influence on the mean flow field. On the other hand, however, the mean flow field 
directly affects the evolution of turbulence through the turbulence dynamical equations. 
This unique feature of the homogeneous turbulence enables us to selectively and efficiently 
study the effect of the mean flow field on the turbulence. 

From turbulence modeling point of view, the homogenous turbulence has been also con- 
sidered as a base for developing advanced closure models, especially the second order 
closure and the spectrum closure (or two-point closure). For the second order closure (or 
second moment closure), the homogeneous turbulence greatly simplifies the turbulence 
equations by excluding the turbulent diffusion term, and thus enables us to concentrate 
on the velocity-pressure gradient correlation tensor and the dissipation rate tensor, which 
are two of the most important terms in the Reynolds-stress equations. In addition, the 
homogeneous turbulence provides us with some nice symmetry properties for turbulence 
correlations which are very helpful in the model development. 

In this paper, we describe several representative second order closure models for the ho- 
mogeneous turbulence, which are also used or extended to the inhomogeneous turbulence. 
Thirteen models! 1 ] - ! 12 ! are included in this study. To explore the behavior of these models 
in a wide range of flows and to identify the ones which represent the state of the art, we have 
made direct comparisons of these models with the data of direct numerical simulations. 
That is, we use the solution of the DNS to calculate the exact values of the terms which 
need to be modeled and then compare them with the values from the respective model 
formulations. This kind of comparison removes the ambiguity and confusion present in a 
global indirect comparison, which compares the solutions of the modeled second moment 
equations. In this indirect comparison, the models of the other terms in addition to the 
model we want to evaluate are also involved. In practice, we often encounter the situations 
in which two wrong things (models) compensate each other to make the solution of the 
specific flow correct, and therefore, we cannot draw a definite conclusion about the certain 
models from this kind of indirect comparisons. 

The DNSs used in this paper are homogeneous shear flows (Rogers et al.! 13 !) and homoge- 
neous irrotational strain and relaxation flows (Lee et al.! 14 !). The DNS data of forty five 
different turbulent flows are used for the direct comparisons of turbulence models. 

Section 2 lists the basic equations and model terms, and then briefly describes the model 
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constraints and the modeling procedure. Section 3 describes various second order closure 
models proposed by different researchers. A brief comment on each model is also given. 
Section 4 shows the detailed direct comparisons between models and DNS data. Finally, 
a few concluding remarks are made in section 5. 

2. Turbulence equations and model terms 

For incompressible turbulent flows with constant density p, the mean equations can be 
written as: 

Ui,i = 0 ( 1 ) 

^r-Ui = -~P,i - (uiu j)j + fcQ + uU tJj (2) 

Ut p 

^0 =-(%),, +7©.ii (3) 

where ( ),» represents the spacial derivative, and ^ ^ Uii ® an< ^ ^ are 

the mean of the velocity, the scalar and the pressure. Ui, $ and p are the corresponding 
fluctuating quantities. #(= pi/p), v and 7 are the buoyancy vector, the kinematic viscosity 
and the thermal diffusivity respectively. To close above equations we need models for the 
second order moments: UiUj , which represent the Reynolds stress and the turbulent 

scalar flux (e.g., heat flux). At the second order closure level, these second moments are 
provided with a set of modeled transport equations. 

The exact equation for the Reynolds-stress tensor UiUj can be written as: 

= Dij + Pij -F Fij + Tij + Ilij — £ij (4) 


where, 


Dij = v(uiUj), kk 

Pij = -I TiUkUj'k - UjWkUi'k 
Fi 


1 j — fiiOuj + (3jdui 

Vij = ~[u l u J u k + \uip6jk + U]p6ik)},k 


n ij - - p(u itJ +u Jfi ) 
tij = 2 UV~^U^k 

Djj and Tij are the viscous diffusion and the turbulent diffusion; P ij and Fij are the shear 
and buoyancy production terms; and finally Ilij and tij represent the pressure-strain rate 
tensor and the dissipation rate tensor respectively. 6ij is the Kronecker’s delta. Obviously, 
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to close above Reynolds-stress equation, we must model the new unknowns T;j, and 
£ij. At the level of the second order closure, these new unknowns are usually modeled 
with algebraic equations in terms of the second moments and the mean quantities (with 
the exception of the trace £kk = 2 e, which is modeled with a transport equation). 

The exact heat flux 0Ui equation is: 

D 

— = Dio -f Pie 4- Fie + Tie + n»0 — (5) 


where, 

D i0 = 7 (uiO tk ),k + v(ui t k0),k 
Pie = ~0UkU i<k -uJnQ,k 
Fe = PjP 

To = —\0uiUk + -p05ik],k 
P 

Ilio = -p0,i 

P 

tie = (1/ + 7)0 

The physical meaning of above terms is similar to the terms in the Reynolds-stress equation. 
If 7 = i/, Die can be written as v(0Ui)kk • Again, in order to close the heat flux equation, 
Tie i Ilifl and must be modeled and they are usually modeled with algebraic equations 

in terms of the second moments and the mean quantities. 

In general, we also need an equation of the temperature variance 0 2 : 

- 2 ^ 0 , * - (0hTk),k - 2 1 0^r k ( 6 ) 

Again, the last two terms in the above equation must be modeled. 

In a general turbulent shear flow with moderate inhomogeneity, the turbulent diffusion 
terms in the second moment equations are usually smaller than the other terms. How- 
ever, the pressure-strain rate and dissipation rate tensors are always among the leading 
terms. Therefore, the performance of modeled equations largely depends on the models of 
pressure-strain rate tensor and dissipation rate tensor. 

For the homogeneous turbulence with constant gradients of the mean velocity and the mean 
temperature, the equations (1-3) will be decoupled from the equations of second moments 
(4), (5) and ( 6 ), and hence the mean flow field of the homogeneous turbulence will not be 
affected by the turbulence. However, the converse is not ture. The mean flow field (Ui 
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and 0) will directly affect the evolution of turbulence through the terms on the right hand 
side of the second moment equations. In addition, the turbulent diffusion terms will not 
be present in the homogeneous turbulence, which enables us to isolate the pressure-strain 
rate tensor and the dissipation rate tensor from the complicated triple correlation tensor 
(turbulent diffusion term). 


In this paper, we will only concentrate on the models of the pressure-strain rate tensor and 
the dissipation rate tensor for the velocity field. However, for the purpose of describing 
the modeling procedure, we will also include the scalar field. 

Realizability: 

For constructing turbulence models, various model constraints have been proposed by dif- 
ferent authors in an attempt to make the model equations as general (or universal) as 
possible. Besides the conventional model constraints (e.g. invariance principle), the most 
recent ones are the realizability principle (Schumann^ 1 ^, Lumley^), the linearity principle 
(Popel 16 J), the rapid distortion theory (Reynoldst 17 !), and the material indifference princi- 
ple (Spezialet 18 !). However, some of the above mentioned principles are not universal. For 
example, the material indifference is not valid for general turbulence in which the fluctu- 
ating velocities are three dimensional. The principle of linearity is also not universal, as it 
holds only for passive scalars. On the other hand, realizability (defined as the requirement 
of the non-negativity of turbulence energy components and Schwarz’ inequality between 
any fluctuating quantities) is the basic physical and mathematical principle that the so- 
lutions of any governing equations should obey. Hence, among all the above mentioned 
model constraints, realizability is the most universal, important principle and is also the 
minimal requirement to prevent the turbulence model equations from producing unphysi- 
cal results. Realizability can be applied to various turbulence quantities, for example, the 
Reynolds stresses v~uj, the scalar fluxes $ul and the triple correlations. For one-point sec- 
ond moment equations, realizability for uiu~j and Oui is most important (for more details, 
see Shih et alJ 19 ), Lumleyf 23 !). 


Modeling procedure: 


A traditional way to treat pressure-strain rate tensor was first proposed by Chout 20 i. It 
starts with the following equation of the fluctuating pressure: 


~ Pijj ~~ u i,j u j,i 

P 


Based on the linearity of p, this equation can be split into three parts: 

-~Pjj = 


(V 

(8) 
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( 9 ) 


J -pf s = -Mi m 

Then the pressure related correlation terms in the second moment equations can be written 
as: 

n« = nSj* + irlf + 14 3) (ii) 

n» = n<J> + + ng> (12) 

where, the first term on the right which is explicitly related to the mean velocity gradient 
is called the rapid term, the second term which is related only to the pure fluctuating 
quantities is commonly called slow term, and the third term which is directly related 
to the buoyancy is called buoyancy term. In the literature, the slow term 7rJ. is often 
combined with the deviatoric part of the dissipation rate tensor (the dissipation is also 
related only to the pure fluctuating quantities). They both drive the turbulence towards 
isotropic state. Therefore Lumley^ 7 ! defined and 

2 

= 7T^ 2) - Sij + -etfij (13) 

- e i0 (14) 

r 

which are called the return-to-isotropy terms in the second moment equations. Here q 2 = 
£ — 2^^^* 

For the rapid and buoyancy terms, one may obtain their exact expressions for homogeneous 
turbulence using the solution of Eq.(8) and Eq.(10): 


= -2U p , q ^- J [(u,(r)u<(r')) iW - + («,(»>,■(»•')),, 


— ^Up^q^Xpjqi *b Xpi q j ) 


nS]’ = -2 J MrWrOL 


2 Uj'kXijk 


n « > = A S / + («(r)»i(r')).«] 

V 

— —0k{Ykji + Ykij) 
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n< ‘« ,=A s/ Wr)tf(r ' )1 ’ i ‘F37i 

V 

= -ftTife (18) 

where the tensors X and Y are the integrals of two-point correlations over the entire 
physical space. To model the rapid and buoyancy terms is now to model these tensors. 
We notice that the exact expressions of X and Y in Eqs.(15-18) do not provide the models 
in terms of the second moments. However, they do give us a hint to model them, because 
these integrals tell us some important properties of X and Y : 


Xpjqi — Xjpqit 

X{jk — Xjik 
Y kji = Y jki 


Yik = Y ki 


^pjqi 


= X 


PJ'Q 


(19) 


Xpjpi = 0, Xikk — o 

Y kjj = 0 (20) 


Xppqi — “UqUi) X{i k fl'U'k 

Y k ki = Ykk=P ( 21 ) 

Eq.(19) is the symmetry condition, Eq.(20) is the incompressibility condition, and Eq.(21) 
is called the normalization condition (Rottat 6 )). Theses equations are very helpful for 
constructing models of X and Y . In fact, they have been used in all conventional turbulence 
models. A model which does not satisfy the above conditions has little hope of success in 
general applications. 

In this paper, we only present the models for the velocity field which will be directly 
compared with the available DNS data. These models are for the rapid term H^ and the 

return-to-isotropy term — e$ij (denoted as ). 

3. Closure models: 

In this section, we will describe five models for the rapid term and eight models for the 
return-to-isotropy term. These models were proposed by different researchers and have 
been indirectly tested in different selected flows. 
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3.1 Models for the rapid term 


Launder, Reece and Rodi (LRR); M 


5% = 0.2 Sij + *Q ± *(b ik S jk + b jk S ik - hijbuSu) 

2 g 2 2* <3 

10 — 7 C 3 ,, n ~ x 

H 22 {bik^ljk 4- bjkQik) 


( 22 ) 


where C 3 = 0.4, and 


s» = \(Ui,j + Uj,i), 

n« = i(Pij - Pj.i) 


This model is linear in the Reynolds-stress. It contains only one model constant C 3 . This 
model satisfies the conventional model constraints Eqs. (19-21). It is the most general form 
at the level of linear dependence on the Reynolds-stress. However, as Lumleyf 7 ! pointed out 
that this model may violate realizability as turbulence approaches two component state. 


Speziale, Sarkar and Gatski (SSG); t 2 l 


n^ 1} 

2 ? 


0.2 -cr c c;p u 

Oij — - 


2g 2 


C 2 

+ -j(h«.Sjk + kjkSn - ^SijhuSu) 


4 - — (bikfljk + bjkftik) 


(23) 


where, 

Cl* = Clbijbij, P = —uiujUij 
Cl = 1.8 Cl = 1.3, C 4 = 1.25, C 5 = 0A 

This model is quasi-linear in the Reynolds-stress, because the coefficients in the first two 
terms are not constant, they depend on the invariant of the Reynolds-stress tensor and the 
production P. This model contains four model constants (Cj*, Cl , C4, C5), therefore one 
may imagine that it will be difficult to correctly calibrate them. In addition, this model 
does not satisfy the basic model constraint Eq.(21). If we impose this constraint, then the 
four coefficients will reduce to only one, and this model will reduce to the LRR model. 
Finally, like the LRR model, the SSG model may also violate realizability. 
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Fu, Launder and Tselepidakis (FLT): ^ 


[\V 2 

4=- = 0 .2Sij + 0.3(bikSjk + bjkSik - -SijbkiSki) 

lq2 o 

13 

H Q~{bikQjk 4" bjk^ik) 

-f 0.2 + bjiSn — 2b kj b li Ski — SbijbkiSki) 

+ o.2 (blciji + b 2 jt n it ) 

+ r[46 nn(bikftjk + bjkSik) 

"f* 126nit^nj(^mfe^njlt ” 1 " 


( 24 ) 


where r = 0.7, b ?- = bikbkj. 

This model is cubic in the Reynolds-stress. The final selected form contains one model 
constant. This model only satisfies a part of realizability condition, that is the two com- 
ponent state of turbulence. However, when scalar field is involved, this model will not 
be able to satisfy the Schwarz’ inequality between velocity and temperature. This part of 
realizability is sometimes called the joint realizability. 

Shih and Lumley (SL): W 


Y=r = 0.2 Sij + 3as(bikSjk + bjkSik — -SijbkiSki) 

+ ^(2 — 7 a 5 )(bikfljk + bjk^ik) (25) 

+ 0.2(6^ Sji + b 2 t Su - 2bkjbuSki — SbijbkiSki ) 

+ 0.2(b 2 [Qji 4- bj/flu) 

where, 

= — (1 + 0.8.F 1 / 2 ), F = 1 + 9bijbjkbki — - bijbij 

This model is quasi-quadratic in the Reynolds-stress, because the model coefficient 05 is 
the function of the invariants of Reynolds-stress tensor. We emphasize that this model 
is obtained from a more general form of the expression than the FLT, and satisfies both 
the two component condition and the Schwarz’ inequality between the velocity and scalar 
fields. In addition, the final form is simpler than the model of FLT. 
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Shih and Mansour (SM); t 5 1 

n* 1 * 2 

— =- = 0 . 2 Sij 4 - 3 a 5 (bikSjk + bjkSik — - SijbkiSki ) 

2 q 2 & 

+ -(2 — 7ccs)(bik£ljk + bjkSik) (26) 

+ 0.2(b^Sji + VjtSu — 2bkjbuSki — SbijbkiSki ) 

+ + tf-ftu) 

where, a 5 = ^{1 + 3.5[1 - (1 - F) 1 / 4 ])}. 

This model has the same form as the SL model. It was derived in a different way and 
contains a different model coefficient as which was calibrated from one of the DNS data 
(Rogers^ 13 !). This model, like the SL model, fully satisfies realizability conditions. 

3.2 Models for the return-to-isotropy term 11^: 


Rotta; ! 6 ! 

n(f = -eCbij (27) 

where, C = 3.0. 

This model is linear in the Reynolds stress, and contains one model constant. It was 
widely used and adopted in the LRR model. We notice that this model will not allow the 
turbulence to reach the state of two component, because when any turbulent component 
reduces to q 2 / 9, the model Eq,(27) will force it to grow. 


Lumley; ! 7 ) 

n'f = -ctfhj + -,(b% + 2/I«y/3)] (28) 

where, 7 = 0 and 

0 = 2 + j exp( — 7.77/v / Fe){72 / \fRe + 80.1 ln[l + 62.4(-/7 + 2.3 III)]} 



This model is quasi-linear in the Reynolds stress, because 7 is set to zero, and (3 is a 
function of the invariants of Reynolds stress tensor. This model is simple, and satisfies 
realizability. 


Sarkar and Speziale (SS); M 

n!]’ = -e[Cife, - 3(C, - ml - (29) 
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where C i = 3.4. 


This is a quadratic model in the Reynolds-stress tensor. It satisfies what they call the weak 
realizability condition. Like the Rotta model Eq.(27), this model will not produce unphys- 
ical results. However, it will not allow the turbulence to approach the two component 
state, which could occur in some situations, for example, in the near-wall turbulence. 

Haworth and Pope (HPW 9 1 

ng> = -£{ 0 , 6 * - cylfc, + + s,j/3)]} (30) 

where C\ = 8.3, C 2 = 14.8. 

Eq.(30) is the slow part of the Haworth and Pope’s model for the situations with no mean 
velocity gradient. This model, like the SS model, will not produce unphysical results, 
however, it will also not allow the turbulence to approach the two component state. 

Choi and Lumlev (CL); I 10 1 


If III > 0, 

where, 


(3 = 2 + 


n{f = -epiy + 7 (»?, + 2 //«y/ 3 )] 

p*F J / 2 


(31.1) 


7 = 


1 + Gx 2 

p*F 1 / 2 G 


1 + g x 2 i 

i = (III/2) 1 / 3 , T? = (-/// 3) 1/2 

X = ~» <? = -x 4 + 0.8 X 6 

V 

P- = «p[-9.2»/ifcV>]{(^“. + 1+L) - [296 - 16.2(X + l) 4 ]//} 


Re 


9ev' 


If III < 0, 


Re'/ 2 Re 
II = —bijbij! 2, III = bijbjkbki/3 

n'J’ = £,.(28) 


(31.2) 


The model coefficients in Eq.(31.1) were obtained using their comprehensive measurements 
of turbulence relaxing from axisymmetric expansion. Both Eq.(30.1) and Eq.(31.2) satisfy 
realizability, however, Eq.(31.1) is valid only for III > 0, because ( is not defined when 
III < 0. 
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Craft and Launder (CfcL): t n l 


= — c,e[ 26y + 4Ci(4? j - 6| 4 Sy/ 3)] - 2riy (32) 

where, 

<7, = 3.1(A 2 A) J / 2 , C*i = 1.2 

g 

A 2 = 4 bijbji, A 3 = 8 bijb jle bki, A = 1 - ~(A 2 - A 3 ) 

O 

This model is tensorially quadratic in the Reynolds stress, and satisfies realizability. 
Yamamoto and Arakawa (YA); t 12 l 

n!j> = -sMy + a 2 (bl - iL<«/3)] (33) 


where 

«! = 2 +P F [q ( b\ k y + | bl k \ s sign(b 3 kk )} 
a 2 = 3 (ai - 2) 

p = -12, q = -0.65, r - 0.4, s = 0.45 

F= + 

The YA model tried to fit the situations with both positive and negative b kk . It also meets 
the requirement of realizability. 

Shih and Mansour (SfcM); f 5 l 

njf = -£{(2.0 + C,F<-)b xj + 1% + (1/3 + 2II)b„ + hl6„\\ (34) 

where, 


Cf = (1/9) exp( — 7.77/\/Re){72/ \/lte + 80.1 ln[l + 62.4(-// + 2.3 III)]} 

7 = To(l — F v ), Re= — 

9ev 


F = 1 + 911 + 3/7/ 


II = -\bnb 


1 


2 III — nbijbjkbki 

( = 17/20, t) = 1/20, 7 o = -2 


This model matches the data of Comte-Bellot and Corrsint 21 ^ and meets the requirement 
that there will be no return to isotropy in the zero Reynolds number limit. This model 
also satisfies realizability. 
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4. Direct comparison between models and DNS data 

The direct numerical simulations of the homogeneous turbulence used for the direct com- 
parisons in this paper were done by Rogers et al^ 3 ^ and Lee et al^ 14 K The Rogers et 
al.’s flows are homogeneous shear flows with different shear rate S and different turbu- 
lent Reynolds number Four of them are marked as C128U, C128V, C128W 

and C128X. They are used for direct comparisons with both rapid and return-to-isotropy 
models. The Lee et al.’s flows are irrotational strain flows and relaxation flows from var- 
ious irrotational strains. The irrotational strain flows include axisymmetric contractions 
(AXK-AXM), axisymmetric expansions (EXO-EXQ) and plane strains (PXA-PXE: the 
cases of suppression in “22” direction, expansion in “33” direction and no strain in “11” 
direction). Eleven of them are used for our direct comparisons with rapid models. The re- 
laxation flows are from axisymmetric contractions (K3R-M5R), axisymmetric expansions 
(03R-Q6R) and plane strains (A2R-H4R). Thirty of them are used for direct comparisons 
with return-to-isotropy models. 

Figures 1-4 show the direct comparisons of five rapid models with the Rogers et al.’s DNS 
data (C128U, C128V, C128W and C128X). The two tensorially linear models (LRRt 1 ' and 
SSGl 2 l) deviate from the DNS data significantly for the II^ and II^ components for all the 
four cases. However they perform reasonably well for the and components except 
for the case of Cl 28V in which the models are off for the component. For most thin 
shear layer flows, the component uju 2 is more important than others. Therefore we may 
expect that the LRR and SSG models would predict thin shear layer flows reasonably well. 
However for flows where all components are important, the above linear models would not 
perform well. We also notice that the SSG model is a quasi-linear model with four model 
coefficients, but according to all test cases, it seems that the SSG model does not show 
any better performance than the LRR model. Apparently, either the model coefficients 
in the SSG model are not chosen properly or the deficiency due to the inconsistency with 
Eq.(21) shows up in these comparisons. On the other hand, three tensorially nonlinear 
models (FLTl 3 ', SLM and SM^ 5 ! ) perform much better than the linear models, expecially 
the SL and SM models, they compare very well with the Rogers et al.’s DNS data for all 
the cases and all the components. 

Figures 5-8 compare seven return-to-isotropy models with the Rogers et al.’s DNS data. 
Two of them (Rottat 6 ! and Lumleyt 7 !) are tensorially linear in the Reynolds stress, others 
are nonlinear. We notice that all the return-to-isotropy models can be written in a basic 
form of Eq.(28) which was proposed by Lumley. In fact, Eq.(28) is the most general form 
provided that is an isotropic function of bij and Re, where (3 and 7 are functions 
of II, III and Re. Therefore, all the return-to-isotropy models are just the variations of 
Eq.(28), depending on the choice of model coefficients (3 and 7 . It is evident from these 
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direct comparisons that the linear Rotta! 6 ! model does not perform very well except for 
the component n^, in which the Rotta! 6 ! model does very well. However this component 
is not very important for most shear flows. The SS model! 8 ! and the YA model! 12 ] are 
nonlinear models; they behave very much like the Rotta model except in the component 
IIjj . On the overall, the SS and YA models perform worse than the Rotta model. The 
C&L model! 11 ! is also a nonlinear model, it behaves better than the Rotta model for all 
the cases except in the component in which the C&L model is the worst one among 
all the models. Apparently, the model coefficients in the above mentioned models are not 
appropriate according to the Rogers et al.’s shear flows. However, surprisingly enough, 
the Lumley’s quasi-linear model! 7 ] performs very well for all the cases. It was well known 
that the Lumley’s model works excellent for the flows in which III < 0, but it would 
not work very well for the cases! 10 ! with III > 0. Here in the Rogers et al.’s flows, 
III > 0, and it still works quite well. We think this is partly due to the low Reynolds 
number behavior of the Lumley’s model and the Rogers et al.’s DNSs are the low Reynolds 
number flows. The CL model! 10 ! j s particularly designed for the cases with III > 0 based 
on their experiments of flows relaxing from the axisymmetric expansion. It also works 
reasonably well in all the Rogers et al.’s flows. Finally, the S&M model! 5 ], which is derived 
from Eq.(28) using realizability and matches the behavior of the low Reynolds number 
turbulence (final period of decay), performs very much like the Lumley’s linear model with 
just a little improvement over it. 

From these direct comparisons, it is clear that the SL! 4 ! and SM! 5 ! nonlinear models perform 
the best among the five rapid models. For the return-to-isotropy models it is also clear 
that the Lumley’s linear model! 7 !, the S&Mt 5 ! model and the CL! 10 ! model are the best 
among the seven models tested here. 

The Lee et al.’s flows are also homogeneous but their characteristics are very different from 
the Rogers et al.’s homogeneous shear flows. Even though these flows do not often occur 
in the nature, it is still interesting to see how do the models of rapid and return-to-isotropy 
terms perform in these critical numerical simulations. We will first look at the comparisons 
of the rapid models and then the return-to-isotropy models. 

Figures 9-11 compare the rapid models with the three axisymmetric contraction flows. 
These figures show that no rapid models, except the SLl 4 l and SM! 5 !, can predict the 
simulation data well. The SSG and LRR models, which are tensorially linear, perform 
the worst. The SSG is even worse than the LRR. The nonlinear model FLT is much 
better than the SSG and the LRR, but still significantly deviates from the DNS data, 
especially in the case AXM, which has high shear rates ( S = 38.2 — 96.5). We notice that 
in the axisymmetric contraction flows, only the SL and SM rapid models show very good 
performance for all the cases and all the components. Figures 12-14 compare the same 
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rapid models with the three axisymmetric expansion flows. This time, no rapid models 
can follow the DNS data well. However, the SM model performs a little better than others, 
and the SSG model is the worst, especially when the shear rate becomes larger. Figures 
15-19 compare the rapid models with the five plane strain flows. The linear LRR model 
does a good job for the component, but bad job for the component. The SSG 
model performs worse than the LRR model. The nonlinear models (FLT, SL and SM) do 
good job for the component, but for the components, they are even worse than 
the linear models. 

Now let us look at the direct comparisons of eight return-to-isotropy models with the Lee 
et al.’s relaxation flows. Figures 20-25 compare the models with the six flows relaxing 
from axisymmetric contractions (K3R-M5R). In these flows, III < 0, the Lumley’s model 
works excellent for all the cases and all the components (the CL model in this case is the 
Lumley’s model). The nonlinear model of S&M works as good as the Lumley’s model. All 
other models deviate from the DNS data significantly, especially when the flows relaxing 
from higher strain rates (e.g. M5R). It is also interesting to note that the simple Rotta 
model works better than the nonlinear models of SS, YA, C&L and HPt 9 l. The HP model 
performs worst among all the models. Figures 26-27 compare the models with the two flows 
relaxing from axisymmetric expansions (03R, 06R), for which III > 0. Figures show that 
the nonlinear models of S&M( 5 1 and CL have better behavior than the other models. This 
time, the model of C&L performs the worst. Finally, Figures 28-49 compare the models 
with the twenty two relaxation flows from plane strains (A2R-H4R). In these flows, the 
Lumley’s linear model and nonlinear models of S&M and CL predict the component 
very well for all the cases. For n^, they are also the best among others. However, this 
time, the SS^ and HPt 9 l models appear to be the best in the component, and the 

YA and C&L models are the worst among others. 

From the above critical comparisons with the Lee et al.’s irrotational strain flows, we see 
that the rapid models of SLM and SMl 5 l work very well in all the axisymmetric contraction 
flows. In the axisymmetric expansion flows and plane strain flows they do not work very 
well. However, on the overall they still perform better than others. The comparisons of 
return-to-isotropy models with relaxation flows show that the Lumley’s linear modell 7j is 
perfect for all the cases with III < 0, so is the model of S&Mt 5 l. For the relaxation flows 
from the axisymmetric expansions and plane strains, the S&M^ and CL^ 10 ! models show 
better performance. 

5. Concluding remarks 

We have made the direct comparisons of five rapid models and eight return-to-isotropy 
models with the direct numerical simulations of Rogers et al.t 13 l and Lee et alJ 14 J Forty 
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five DNS flows are used for the direct comparisons. We notice that the Reynolds number 
in all these simulations is low, and therefore, they may not represent the real turbulence 
in the nature. However, the model terms concerned here are mainly pressure related 
correlations. Eq.(7) indicates that the fluctuating pressure is not directly related to the 
viscosity, hence the pressure related correlation terms may not be directly affected by the 
Reynolds number, especially the rapid term. The return-to-isotropy term which includes 
the deviatoric part of dissipation rate tensor may have some dependence on the Reynolds 
number. According to the above consideration, we think that the direct comparisons with 
the low-Reynolds DNS data are legitimate, although we should keep in mind the possible 
low- Reynolds number effect of the DNS data. 

We have directly compared five rapid models with fifteen DNS flows: four of the Rogers 
et al.’s shear flows, eleven of the Lee et al.’s irrotational strain flows (axisymmetric con- 
traction, axisymmetric expansion and plane strain). Comparing the performance of the 
LRR and SSG models, which are tensorially linear in the Reynolds stress, we conclude 
that the SSG model gives very little improvement over the LRR model. In fact in many 
cases, it is worse than the LRR model. The reason is not very clear. However, we notice 
that the SSG model does not satisfy the normalization condition of Eq.(7) which may be a 
cause for its poor behavior. If we impose the constraint of Eq.(7) on the SSG model, then 
it will exactly reduce to the LRR model. In fact it can be shown that the most general 
form of the rapid model, which is tensorially linear in the Reynolds stress, is the LRR 
model. Therefore, in general, the treatment used in the SSG model would hardly give any 
improvement over the LRR model. A natural way to improve the model is to use a more 
general nonlinear form and more general model constraints. A typical example is the SLW 
model. It starts with the most general form, using full realizability constraints together 
with the other conventional constraints given by Eq.(5)-Eq.(7). The result is a well be- 
haved model. Indeed, from the direct comparisons with the DNS data, the SLW model and 
its variation form of SM! 5 ! model give the best performance in most of the cases. As to the 
FLT( 3 1 model, it is also a nonlinear model. It starts with a tensorially cubic dependance 
on the Reynolds stress with constant coefficients (in general, these coefficients should not 
be restricted to constants). In addition, the two component conditions of turbulence have 
been imposed. However, the FLT model ignores the Schwarz’ inequality. Its final form 
contains two undetermined model constants, but one of them is set to zero. The perfor- 
mance of the FLT model, from the direct comparisons with the DNS data, is on the overall 
better than the linear models, but does not compare with the performance of the SL and 
SM models. So from these direct comparisons of the rapid models, we conclude that the 
SLW model and its variation form SM^l are clearly the best. Having said this we notice 
that, as Reynolds! 22 ! pointed out, any of the above rapid models will not show any effect 
of the rotation on the invariants (//, III) of the anisotropy tensor b{j. This is clearly a 
theoretical deficiency of the current rapid models. A further investigation is needed to find 
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that how serious this deficiency will be in the practice. 


We have directly compared eight return-to-isotropy models with thirty four DNS flows: 
four shear flows and thirty relaxation flows from axisymmetric contraction, axisymmetric 
expansion and plane strain. As was discussed earlier, all the return-to-isotropy models are 
the variation of Eq.(28) derived by Lumleyf 7 !. Therefore the differences in the models are 
due to the different choices of the model coefficients. Two linear models are the Rotta! 6 ' 
and LumleyM (which is quasi-linear in The Lumley’s model satisfies realizability, 

matches the data of Comte-Bellot and Corrsin^^ and the limit of the final period of 
decaying turbulence. It performs perfectly when III < 0. It also compares well with the 
DNS data in which III > 0. The Rotta’s model does not compare with the performance 
of the Lumley’s model. In fact, the nonlinear models of SS, YA, HP and C&L also do not 
compare with the performance of the Lumley’s model. Apparently the model coefficients 
chosen in these models are not appropriate. The CL^ 10 ! model is designed for the flows with 
III > 0 and is based on their experiments of the relaxing turbulence. It does work better 
than the Lumley’s model when III > 0. Finally, the S&Mt 5 ! model is a nonlinear model, it 
works just like the Lumley’s model when III < 0. When III > 0, it shows an improvement 
over the Lumley’s model according to the DNS data. So from these direct comparisons of 
the return-to-isotropy models, we conclude that the combination of the Lumley’s model 
and the Choi’s model, that is the CLt 10 l model, will give the best performance. The S&M( 5 1 
model seems as good as the CL model according to these comparisons. Having said this, we 
notice that the existing return-to-isotropy models do not follow the relaxation flows from 
expansion and plane strain very well. Therefore there is still a need to further investigate 
and improve the return-to-isotropy models. 
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Figure 1. Direct comparison of the rapid models with the DNS data of the 
homogeneous shear flow C128U (Rogers et alJ 13 !). 
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Figure 2. Direct comparison of the rapid models with the DNS data of the 
homogeneous shear flow C128V (Rogers et alJ 13 l). 
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Figure 3. Direct comparison of the 
homogeneous shear flow C 


,pid models with the DNS data of the 


8W (Rogers et alJ 13 l). 
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Figure 4. Direct comparison of the rapid models with the DNS data of the 
homogeneous shear flow C128X (Rogers et atJ 13 l). 
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Figure 6. Direct comparison of the return- to-isotropy models with the DNS data 
of the homogeneous shear flow C128V (Rogers et alJ 13 ^). 
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Figure 7. Direct comparison of the return-to-isotropy models with the DNS data 
of the homogeneous shear flow C128W (Rogers et alJ 13 J). 
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Figure 8. Direct comparison of the return-to-isotropy models with the DNS data 
of the homogeneous shear flow C128X (Rogers et atJ 13 l). 
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Figure 9. Direct comparison of the rapid models with the DNS data of the 
axisymmetric contraction AXK (Lee et alJ 14 J). 


Figure 10. Direct comparison of the rapid models with the DNS data of the 
axisymmetric contraction AXL (Lee et alJ 14 l). 
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Figure 11. Direct comparison of the rapid models with the DNS data of the 
axisymmetric contraction AXM (Lee et alJ 14 ^). 


Figure 12. Direct comparison of the rapid models with the DNS data of the 
axisymmetric expansion EXO (Lee et alJ 14 !). 
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Figure 13. Direct comparison of the rapid models with the DNS data of the 
axisymmetric expansion EXP (Lee et alJ 14 l). 




Figure 14. Direct comparison of the rapid models with the DNS data of the 
axisymmetric expansion EXQ (Lee et alJ 14 l). 
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Figure 20. Direct comparison of the return-to-isotropy models with the DNS dat 
of relaxation from axisymm. contraction K3R (Lee et al.( 14 l). 
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Figure 21. Direct comparison of the return- to-isotropy models with the DNS data 
of relaxation from axisymm. contraction K6R (Lee et alJ 14 l). 
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Figure 22. Direct comparison of the return-to-isotropy models with the DNS data 
of relaxation from axisymm. contraction L3R (Lee et alJ 14 l). 




Figure 23. Direct comparison of the return-to-isotropy models with the DNS data 
of relaxation from axisymm. contraction L6R (Lee et alJ 14 ^). 
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Figure 26. Direct comparison of the return- to-isotropy models with the DNS data 
of relaxation from axisymm. expansion 03R (Lee et alJ 14 l ). 


Figure 27. Direct comparison of the return-to-isotropy models with the DNS data 
of relaxation from axisymm. expansion 06R (Lee et alJ 14 !). 
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Figure 28. Direct comparison of the return-to-isotropy models with the DNS data 
of the relaxation from the plane strain A2R (Lee et alJ 14 !). 
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Figure 39. Direct comparison of the 
of the relaxation from t’ 
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Figure 47. Direct comparison of the return-to-isotropy models with the DNS data 
of the relaxation from the plane strain H2R (Lee et alJ 14 i). 
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